Vacancy-Impurity Nanoclusters in Solid Solutions 3 He — 4 He 
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The existence of vacancy-impurity clusters due to quantum properties of vacancies in phase 
separated solid solutions of 4 He in 3 He is analyzed and discussed. Additional mechanism called 
vacancy assisted nucleation is proposed. According to this assumption the vacancy-impurity clusters 
should have b.c.c. structure. 

PACS numbers: 67.80 



I. 



INTRODUCTION 



Vacancy-impurity clusters (VIC) were predicted in 
1978 by this author (see also Refs. 0, 0) as one of 
a series structural effects that could appear due to the 
quantum nature of vacancies and impurities in quantum 
crystals. This effect was confirmed in 2001 Q and fur- 
ther investigated 0, 0, in the famous experiments of 
the Kharkov group, using new methods based on a pre- 
cise measuring of the pressure variations in phase separa- 
tion of solid solutions 3 He — 4 He. The effect is based on 
the fact that the derealization of a vacancy in a quan- 
tum solid is attended with a lowering of the energy by an 
amount of the half vacancion band width. This requires, 
however, a good periodicity of the surrounding area. In 
helium solid solutions the periodicity can be violated by 
isotope impurities. Hence, the vacancy could i) push out 
the impurities, ii) attract them and create a cluster of 
impurity atoms, or iii) rearrange host atoms and impu- 
rities in a periodic lattice. We shaw that all three pos- 
sibilities can be realized in the helium solutions. Other 
possibilities can be found in Refs. Ganshin et al. 
IE observed clusters of 4 He impurities around 
vacancies in dilute solutions of 4 He in 3 He. They called 
them Andreev-Pushkarov nanoclusters having probably 
in mind some analogy between the structural effect in 
solid solutions and the magnetic vacancy in pure 3 He 
proposed by Andreev but unfortunately not observed 
yet. 

The aim of this article is to consider some peculiarities 
of the vacancy- impurity clusters and turn attention to an 
additional mechanism in their creation. 



II. THE RID EFFECT 

Let us consider an isolated vacancy in a dilute solution. 
Let Sq be the formation energy of a localized vacancy. 
The derealization lowers this energy by an amount of 
A/2 where A is the vacancion energy band width. Since 
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this is only possible in a periodic lattice, the impurities 
should be remouved from a region of radius R (measured 
in this work in interatomic distances a — (Im/AT^) 1 / 3 ). 
Due to the Heisenberg uncertainty principle the lowest 
energy level in such a potential region is approximately 



e = % h /MR = ir z A/R 2 



(1) 



where we used the relation between the vacancion effec- 
tive mass M and the tunnelling amplitude A for cubic 
crystals M = h 2 /(Aa 2 ) with the interatomic distance 
a = 1. In the same notation the band width is A = zA, 
z being the number of the nearest neighbours (for a h.c.p. 
lattice M = h 2 / (2Aa 2 ) , and the energyband consists of 
two partially overlapping parts |2, EJ The ridding 
changes the entropy. Hence, the change of the free en- 
ergy is of the form 



F = E -TS = e 
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where S — xh\(e/x) is the entropy per unit lattice site, 
and x is the impurity concentration. Expression has 
a minimum at [l| 
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and the volume of the rid region is 



3 V 2T5 / 



(3) 
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This result is restricted by the requirement that the free 
energy be not smaller than £o — A/2, as well as by the 
condition R > a. The latter condition is stronger, and 
reads TS < A. In addition, the energy of the first quan- 
tum level JTJ should be less than A/2. 

Of course, expressions © - O are only estimations 
showing the order of magnitude. These formulae do not 
take into account the change of the impurity concentra- 
tion in the solution. The important result is that the size 
of the rid region is determined by the interplay between 
the amplitude A (the vacancy exchange integral, not the 
band width) and the entropy term and does not depend 
on the "won" energy due to the delocalization, A/2, and 
the impuriton characteristics. The quantum nature of 
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the impuritons shows itself in their mobility. The large 
value of the entropy term is the reason the zone motion 
of defects in 3 He to be suppressed. 

Analogous formula has been obtained by Andreev Q 
for the number of aligned spins N s around a vacancy 
in solid 3 He. In that case S = nln2 = In 2 and the 
expression (J2J takes the form 
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(5) 



This equation with T = 150 mK and A = IK yields 
R 1.5 and N s « 14. Note that with this value 
of R the derealization energy ir 2 A/R 2 w 4.5 K turns 
out to be larger than the half band width A/2 in 3 He, 
i.e. the first quantum level does not lie in the well. 
The minimum value of the free energy turns out to be 
higher than the formation energy of a localaized vacancy: 
F mm = e - A/2+11. 85 A 3 ^ 5 T 2 / 5 -> e - A/2+5. 55 > e . 
Therefore, the effect considered should take place eventu- 
ally in materials with very large bandwidths, but hardly 
in 3 He. 

Let us turn back to the solid solutions 4 He — 3 He. The 
approximation above can be improved 0, for the more 
realistic case of vacancies in a solutiowith a given small 
concentration x v — n v /N and a concentration of impu- 
rities x — n.i/N. As a result of ridding, the impurity 
concentration increases to the value 

(6) 



N - Nx v V 1 - x v V 

where V is the rid area around a vacancy. The problem 
can be formulated as looking for the minimal work for 
changing the concentration in the solution x — ► x' by 
removing x v VN particles of the solvent. 

The Gibbs potential for a solution with concentration 
x is of the form fijl 



$i = iV> + NxT In - + Nxip 
e 



(7) 



where fj,o is the chemical potential of a pure solvent, and 
^ = ^(T,P). 

After ridding the impurity concentration becomes x', 
since the same number impurities is distributed over 
N' = N - x v VN = N(l - x v V) = Nx'/x lattice sites. 
Therefore, the Gibbs potential $2 of the solution has 
the same form J7J) with x and N replaced by x' and 
N' = Nx'/x, while the potential of the pure solvent re- 
moved is $0 = XyVN/iQ- Their sum is 



2 = $ 2 + $ = Nfi + NxTln— +Nxip (8) 

e 



The difference <J> 2 - $1 is 



The change in the energy due to the presence of Nx v 
vacancions in the system is 



A ir 2 A 
AE = Nx v [e -- + 1 p 



(11) 



and finally the free energy per particle takes the form 

F = x v (e - y + ^fj - xT\n (l - x v ^ttR^ (12) 
The free energy F has its minimum at 0, Q 
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In i|13fl the concentration x' stands instead of the entropy 
S = x\n(e/x). The minimal value of the free energy is 
then 



min ^0 
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where n x = xV is the number of impurities rid by one 
vacancy. The change in the free energy affects the equi- 
librium vacancy concentration: 



x v = exp ■ 



£0 - A/2 IOtt (kA\ 
f 3~ \ Tr) 



3/5 
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Hence, the equilibrium vacancy concentration decreases 
with increasing number of defects removed. 

Let us evaluate the volumes Vo and V for the typical 
values A = IK, T = 0.1K, and x = 1 %. Then one 
has by means of eqs. (gj and (SJ V = 123, V = 346. 
For T = 0.2 K and x = 2 % the same expressions yield 
Vo = 60 and V = 150 (we did not take into account an 
eventual dependence of x' on T) . The dependence of the 
number of atoms in a cluster on temperature (at constant 
x') corresponding to formulae @ and (|13fl are shown in 
Fig. 1. In fact, the number of the removed impurities by 
a vacancy n x = xV ~ 2 ~ 5 and the relative change of 
the concentration Sx/x is of the order of several percents. 
However, this can lead to a noticeable decreasing of the 
vacancy equilibrium concentration by more than 2 orders 
of magnitude (see eq. JTBJ). It is worth noting that the 
region considered is not a cluster. It consists of solvent 
atoms and has no boundaries, surface energy etc. 



A$ = NTx In — 

x 



Hence, the variation of the entropy equals 



AS = Nx In ■ 



(9) 



(10) 



III. DIFFUSION 

The ridding effect may be considered as an effective re- 
pulsion which prevents impurities and vacancies be closer 
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FIG. 1: The number of atoms Vo as versus temperature at 
different concentrations, solid line: x' = 1 % eq. HK-^ . dashed 
line: x' — 2% eq. long-dashed line: x' = 2% eq. Q). 



to each other. This means that an effective cross-section 
of impuriton scattering on vacancies can be introduced. 
It can be evaluated by means of the radius R and is of 
the order of a ~ nR 2 ~ 7r(7rA/2Ta;) 2 / 5 . In order to avoid 
any misunderstanding we would like to note that this is 
not the case of 4 He impurities in 3 He, because their zone 
motion is embarrassed by the chaotic orientation of the 
nuclear spins of the lattice atoms. It could take place in 
a strong orienting magnetic field. The values of a for a 
hypothetic quasiparticlc in a medium with helium char- 
acteristics are a ~ 60 a 2 (at T = 100 mK and x = 1 %). 
The temperature and concentration dependencies of the 
diffusion coefficient in this case are 



D = D - 



T 2/5 
,3/5 



(16) 



We considered till now concentrations less than the 
saturation concentration x s . If x' > x s then a new phase 
of impurities appears. The effect of vacancies on this 
phase is considered in the next section. 



IV. VACANCY-IMPURITY CLUSTERS 



The chemical potential n\ of the pure phase corresponds 
to some saturation concentration, x s (otherwise, the 
phase should solve). Hence, fx® = Tlnx s + ip, and 



-SnT In ■ 



-NT5x\n — . (19) 



We suppose that all 5n atoms are confined to n v vacan- 
cies. Then, repeating the procedure in Sec. II yields 
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(20) 



with the restrictions T\ \n(x/x s )\ < A and x ^ x s . If 
the concentration x equals the saturation (equilibrium) 
concentration then the entropy term vanishes and the va- 
cancy is fully delocalized. This effect should be seen in 
a rapid increase of the vacancy mobility inside the pure 
phase. Such an observation was reported in [Tl| where 
fast nonphonon excitations were found in solid 4 He with 
a small concentration of 3 He. Let us mention that x 
is the concentration in the matrix after the clusters are 
formed. The values of x and x s should be taken from 
the experiment and the phase diagram. A quite sophis- 
ticated analysis using known data for helium solutions 
was recently made in Ref. Q • The authors found for the 
cluster size R « 2.1 -=- 4.4 a. A small correction to these 
values could appear if the band width A in the expres- 
sion for the localization energy ir 2 A/R 2 is replaced by A 
(sec cq. 0). Otherwise, the first quantum level e turns 
out to lie out of the potential well with a depth A/2 
(tt 2 A/R 2 > A/2 for all values of R obtained). 

For a rough evaluation one can simplify (|20|) . If the 
initial concentration is xq, then V-n s /Nr* 
x s — x v V where n s is the number of atoms in the pure 
phase and we neglect the small number of solvent atoms 
solved in the pure phase. The initial concentration drops 
out as it should be in a case of phase separated mixture 
in equilibrium. Then (|20|) can be rewritten in the form 



Let us consider a solid solution of impurity atoms (say 
4 He) with a concentration x = n/N where N is the num- 
ber of solvent atoms (say 3 He). If x < x s , Gibbs potential 
is of the form J7J) 



$i = Nu Q + nT\n 



Ne 



nip 



(17) 



where /i refers to the pure solvent and ip = ip(P,T). 
The chemical potentials ^3 and U4 that correspond to 
the solvent atoms and impurities in the solution are 



U3 = fig — Tx, Hi — Tlnx + tp 



(18) 



Let us move Sn impurities from the solution to their 
"pure" concentrated phase. The change of the thermo- 
dynamic potential is 



A$ 



-Sn(T In x + ijj) + Sri/i® 
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and therefore, 
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(21) 
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(22) 



For T = 150 mK, x v = 10~ 5 , x s = 5 x 10~ 3 one has R w 
2.5 and V « 60 (x v V/x s = 0.12 <C 1). Unfortunately this 
approximation requires x v V x s and is not applicable 
to the most interesting low temperature case where x s 
can be smaller than x v V. 

It is worth noting, that VIC can appear both inside an 
existing phase and as a separate inclusion. This makes 
difficult to evaluate the role of the surface energy. 
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V. LOCAL PHASE STRATIFICATION EFFECT, 
VACANCY ASSISTED NUCLEATION 

This effect was predicted in the same work (see 
Refs. H as weii) and is another indication for the 
quantum nature of vacancies. It appears when the va- 
cancion band widths in two phases of one and the same 
substance (e.g. h.c.p. and b.c.c. 4 He) differ significantly. 
Let the h.c.p. phase be the stable one (near the tran- 
sition to b.c.c). If a vacancy appears, the energy in- 
creases by eh — (subscript h corresponds to h.c.p.). 
If it were in a b.c.c. phase, the energy change would be 
efe — Ab/2. Since near the phase separation the thermo- 
dynamic potentials of the pure phases are almost equal, 
the more favorable state is determined by the difference 
SE = e h - A h /2 - (e b - A b /2) = Se + (A b - A h )/2. For 
large enough A b the difference SE > independently of 
the sign of Se and the b.c.c. phase with a vacancy has 
lower energy. Of course, the bottom of the band may 
become lower due to the large difference of the vacancy 
formation energies, but this seems not to be the case 
of helium crystals. The vacancy will rearrange the 4 He 
atoms from h.c.p. to b.c.c. hence creating a nuclei of a 
new phase. The size R of the nuclei can we evaluated 
using the procedure applied above. The difference of the 
free energy is 

5F = 5e+ ^ + TSS-nR 3 (23) 
R A 3 

where A is the vacancion exchange integral in the b.c.c. 
crystal. Having in mind that TSS = q is the latent heat, 
one obtains the nuclei radius 

Let us note, that the surface energy in this case is smaller 
than in the case when the matrix and the nuclei have 
different structure (b.c.c. and h.c.p.). We were not able 
to find reliable data for the latent heat at the temperature 
region considered. If one takes AS ~ 1CP 2 extrapolating 
Pj then q ~ 1CT 3 and R n w 4.4, V n « 350. These 
are probably the clusters observed in the experiments of 
Kharkov group. After the inclusion of 4 He with a vacancy 
appears, a local phase transition to the more favorable 
b.c.c. phase is performed. This makes the VIC even more 
stable. We suppose, that this effect could contribute to 
the curious observation in 0] that the lattice constants 
in both phases of a phase separated solution increase after 
heating. 

Experimental evidences to support our hypothesis 
should be looked for not only in a phase separated so- 
lutions. As we showed in [HUQ inclusions of b.c.c. 4 He 



should manifest themselves in a hysteresis at b.c.c. -h.c.p. 
transitions of pure 4 Hc. In this case one can use the ex- 
perimental data for the transition line h.c.p. -b.c.c. and a 
nuclei contains ~ 200 -j- 250 atoms. Phase nuclei of the 
same kind and of approximately the same size should ap- 
pear in the liquid helium near the liquid-solid transition 
line. 

The theory of VIC is instructive from another point 
of view as well. At the first glance, the small exchange 
integral in spin-spin interaction of 3 He atoms compared 
to the vacancion bandwidth should lead to their easy 
orientation. It does not happen due to the large entropy 
factor. So, vacancy prefers to rearrange atoms (not spins) 
organizing a local phase of impurities and to delocalize in 
this phase. Its choice is the minimal change of entropy. 



VI. CONCLUSION 

In this work three effects due to the quantum nature of 
vacancies in solid helium were considered relative to the 
vacancy-impurity complexes observed recently in phase 
separated solid solutions of 4 He in 3 He. The size of a 
cluster depends on the ratio of the vacancy exchange in- 
tegral A and the entropy term TAS, and does not depend 
directly on the vacancion band width. The experimen- 
tal results cannot give, therefore, information about the 
band width. Nevertheless, it is essential for the phase 
transition. It is supposed that VIC in rear 4 He - 3 He 
solutions undergo two transitions: first complexes with 
h.c.p. structure appear and then they transfer to b.c.c. 
clusters. The transition is controlled by the positions of 
the bottoms of the vacancy bands, and by the minimum 
entropy change. The evaluation of the clusters size is in 
good agreement with the experimental observations. It 
is shown, that near the separation line the size of the va- 
cancy cluster rapidly increases, the vacancy becomes de- 
localized and mobil inside the 4 He- phase. It is supposed 
that this could be the fast nonphonon mode observed in 
[Tl|. The effect of vacancy assisted local phase stratifi- 
cation should be observable near the b.c.c. -h.c.p. as well 
as b.c.c. -liquid transition. 
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